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Abstract 

In this article, I calculate the contributions of the nuclear matter induced condensates 
up to dimension 5, take into account the next-to-leading order contributions of the nuclear 
matter induced quark condensate, study the properties of the scalar, pseudoscalar, vector 
and axialvector heavy mesons in the nuclear matter with the QCD sum rules in a systematic 
way, and obtain the shifts of the masses and decay constants. Furthermore, I study the heavy- 
meson-nucleon scattering lengths as a byproduct, and obtain the conclusion qualitatively about 
the possible existence of heavy-meson-nucleon bound states. 
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1 Introduction 

The suppression of J/ip production in relativistic heavy ion collisions is considered as an important 
signature to identify the quark-gluon plasma [Tj. The dissociation of J/ip in the quark-gluon plasma 
due to color screening can result in a reduction of its production. The interpretation of suppression 
requires the detailed knowledge of the expected suppression due to the J/ip dissociation in the 
hadronic environment. The in-medium hadron properties can affect the productions of the open- 
charmed mesons and the J/ip in the relativistic heavy ion collisions, the higher charmonium states 
are considered as the major source of the J/ip [5j. For example, the higher charmonium states 
can decay to the DD, D*D* pairs instead of decaying to the lowest state J/ip in case of the mass 
reductions of the D , D *, D , D* mesons are large enough. We have to disentangle the color screening 
versus the recombination of off-diagonal cc (or bb ) pairs in the hot dense medium versus cold nuclear 
matter effects, such as nuclear absorption, shadowing and anti-shadowing, so as to draw a definite 
conclusion on appearance of the quark-gluon plasma mm- The upcoming FAIR (Facility for 
Antiproton and Ion Research) project at GSI (Institute for Heavy Ion Research) in Darmstadt 
(Germany) provides the opportunity to study the in-medium properties of the charmoniums or 
charmed hadrons for the first time. The CBM (Compressed Baryonic Matter) collaboration intends 
to study the properties of the hadrons in the nuclear matter [5], while the PANDA (anti-Proton 
Annihilation at Darmstadt) collaboration will focus on the charm spectroscopy, and mass and 
width modifications of the charmed hadrons in the nuclear matter (jj . However, the in-medium 
mass modifications are not easy to access experimentally despite the interesting physics involved, 
and they require more detailed theoretical studies. On the other hand, the bottomonium states 
are also sensitive to the color screening, the T suppression in high energy heavy ion collisions can 
also be taken as a signature to identify the quark-gluon plasma [7]. The suppressions on the T 
production in ultra-relativistic heavy ion collisions will be studied in details at the Relativistic 
Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC). 

Extensive theoretical and experimental studies are required to explore the hadron properties 
in nuclear matter. The connection between the condensates and the nuclear density dependence 
of the in-medium hadron masses is not straightforward. The QCD sum rules provides a powerful 
theoretical tool in studying the in-medium hadronic properties ED 121, and has been applied exten¬ 
sively to study the light-flavor hadrons and charmonium states in the nuclear matter nni mi ra- 
The works on the heavy mesons and heavy baryons are few, only the D , B , Do, Do, D*, B*, D±, 
Bi, A q, Eg, E qq and Qqq are studied with the QCD sum rules 1131 HU. G, SI HZ1 HU- The 
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heavy mesons (heavy baryons) contain a heavy quark and a light quark (two light quarks), the 
existence of a light quark (two light quarks) in the heavy mesons (heavy baryons) leads to large 
difference between the mass-shifts of the heavy mesons (heavy baryons) and heavy quarkonia in 
the nuclear matter. The former have large contributions from the light-quark condensates in the 
nuclear matter and the modifications of the masses originate mainly from the modifications of 
the quark condensates, while the latter are dominated by the gluon condensates, and the mass 
modifications are mild 

In previous works nn ei Ha mum, the properties of the heavy mesons in the nuclear matter 
are studied with the QCD sum rules by taking the leading order approximation for the contri¬ 
butions of the quark condensates. In this article, I take into account the next-to-leading order 
contributions of the quark condensates, and study the properties of the scalar, pseudoscalar, vec¬ 
tor and axialvector heavy mesons in the nuclear matter with the QCD sum rules in a systematic 
way, and make predictions for the modifications of the masses and decay constants of the heavy 
mesons in the nuclear matter. Measuring those effects is a long term physics goal based on further 
theoretical studies on the reaction dynamics and on the exploration of the experimental ability to 
identify more complicated processes mm- Furthermore, I study the heavy-meson-nucleon scat¬ 
tering lengths as a byproduct. From the negative or positive sign of the scattering lengths, I can 
obtain the conclusion qualitatively that the interactions are attractive or repulsive, which favor or 
disfavor the formations of the heavy-meson-nucleon bound states. For example, the £ c (2800) and 
A c (2940) can be assigned to be the S'-wave DN state with J p = i and the S-wave D*N state 
with J p = | respectively based on the QCD sum rules m- 

The article is arranged as follows: I study in-medium properties of the heavy mesons with the 
QCD sum rules in Sec.2; in Sec.3, I present the numerical results and discussions; and Sec.4 is 
reserved for my conclusions. 


2 The properties of the heavy mesons in the nuclear matter 
with QCD sum rules 

I study the scalar, pseudoscalar, vector and axialvector heavy mesons in the nuclear matter with 
the two-point correlation functions 11(g) and II At!y (g), respectively. In the Fermi gas approximation 
for the nuclear matter, I divide the 11(g) and II M „(g) into the vacuum part II 0 (g) and n° y (g) and 
the static one-nucleon part IlAr(g) and II^(g), and expand the Hn{q) and 11))), (g) up to the order 
O(pn) at relatively low nuclear density pTTl fT5j, 


n(g) = iJd*xe* x (T{J(x)j'(0)}) PN 

= n 0 (g) + -^T iv (g), 

2m n 

IWg) = ijd l xe**(T{J li (x)Jt(0)}) PN 

= + « 

where the p jv is the density of the nuclear matter, and the forward scattering amplitudes T)v(g) 
and T/)f, (g) are defined as 

T N (u>,q) = i J d 4 ‘xe ig ' x {N{p)\T , 

T”(u,q) = i j d A xe^^N(p)\T{j^{x)Jt{{))}\N(p)), (2) 
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where the J{x) and J^,{x) denote the isospin averaged currents r)(x), 775(2), 77^(2) and 775^( x ), 
respectively, 


77(2) 

t c(x)q(x) + q(x)c{x) 

= V'(x)= 2 


V 5 {x) 

t 0(2)7759(2) + 7(2)7750(2) 

= %\ x ) = 2 


V»{x) 

t , ^ c(x)l^q{x) + q( 2)7 M c(2) 

= %(x) = 2 


V 5 »{x) 

„t 7~\ 5(2)7^759(2) + 9(2)7^750(2) 

'/s^W — 2 ’ 

( 3 ) 


which interpolate the scalar, pseudoscalar, vector and axialvector mesons D 0 , D , D* and 
respectively. I choose the isospin averaged currents since the Do, D , D* and D\ mesons are 
produced in pairs in the antiproton-nucleon annihilation processes. The 9 denotes the u or d 
quark, the q M = ( 07 , 9 ) is the four-momentum carried by the currents J(x ) and J^{x), the \N(p)) 
denotes the isospin and spin averaged static nucleon state with the four-momentum p = (mw, 0 ), 
and (N(p)\N(p')) = (2ir) 3 2p 0 5 3 (p - p’) [15] , 

I can decompose the correlation functions T^ v (oj,q) as 

T* v {u,q) = T N (u},q)(^-g^ v + ^^j +T^(uj,q)q^q„ 

+T^(w, q ) (q^Uv + q v u^) + T^(w, q )u^u v , (4) 

according to Lorentz covariance, where the T/v(w,q) denotes the contributions of the vector and 
axialvector charmed mesons, and the T^' 1 ^ 2 (uj, q ) are irrelevant in the present analysis. 

In the limit q —>• 0 , the forward scattering amplitude T/v(w, q ) can be related to the DN (Do-ZV, 
D*N and D\N ) scattering T-matrix, 

Td/Dq/d*/d 1 n{vid/Do/d*/D u®) = 87t(toat + rn D / Do / D */ Dl )a D / Do / D */ Dl , (5) 

where the cld/Do/d*/d x are the D/Dq/D*/D\N scattering lengths. I can parameterize the phe¬ 
nomenological spectral densities p(u> , 0) with three unknown parameters a , b and c near the pole 
positions of the charmed mesons D , Do, D* and D\ according to Ref.[I3l. 


p(w, 0 ) =-Im 

7 T 


Td/d 0 n{u, 0) 


lb/D 0 m D/D 0 


+ 


_( u2 - m D/D 0 + ie ) 

= « ^ 2 6 ( w2 - m2 D/D 0 ) +b6(u} 2 - m 2 D/D ^j + cS (w 2 - s 0 ) , 
for the pseudoscalar and scalar currents 775(2) and 77(2), 


( 6 ) 


p(w, 0) =-Im 

7 T 


Td*/d 1 n{u, 0 ) 


( 


u 2 -m‘ 2 D „/ Dl 


■ l£ 


TD*/D 1 ni D* /Dx 


= a J^ S ( u) 2 - m D-/Di) + b $(u 2 -m 2 D , /Dl ) +cS(u 2 -s 0 ) 
for the vector and axialvector currents 77 ^( 2 ) and 775 ^( 2 ). 


( 7 ) 
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Now the hadronic correlation functions n(w, 0) and 0) at the phenomenological side can 

be written as 


n(w,o) = 


( Id/Do + fifD/D 0 ) (®D/D 0 + ^ m D/Do) 


1 


f 2 777 4 

J D/D 0 ln D/D 0 


1 


(mD/Do + Sm D/Do ) - w 2 


mi 


m D/ Do ~ w2 


+ • 


Pn 

2 toat 


+ 


0 m D/D 0 - W 


2 V m 2 D/D 0 U 2 


( 8 ) 


n M „(w,0) — (/d*/Di + <5/d*/£>i) (jnD*/D! + 5m D */ Dl ) 


<7/7 91. 

<?/ 7 !/ + 2 


- f D* /Di ^D* /D\ 2 


m D*/D,-“ 


Pn 

2 toat 


( 


u 


m2 D^/D 1 - Uj2 


9/7 917 

2 I 9/717 + 2 


(m D */ Dl +Sm D . /Dl ) - w 2 


1 2 m D*/.Di w2 


, 9/7 9j- 

9/717 + 2 


(9) 


In Eqs.(6-7), the first term denotes the double-pole term, and corresponds to the on-shell effect 
of the T-matrix, 


= - 8 tt (tow + m. D jD 0 )a D /D 0 


D/D 0 rn D/D 0 

to 2 


( 10 ) 


for the currents 775 ( 2 ;) and 77 ( 2 ;) and 


= -87r(mAr + m D , /Dl )a D * /Dl fl, /Di m 2 D , /Di , 


( 11 ) 


for the currents 77,7(2;) and 775,7(2;); the second term denotes the single-pole term, and corresponds 
to the off-shell effect of the T-matrix; the third term denotes the continuum term or the remaining 
effects, where the so is the continuum threshold parameter. In general, the continuum contribu¬ 
tions are approximated by pqcd{w, 0 )9{ur — s 0 ), where the pqcd{w, 0) are the perturbative QCD 
spectral densities, and 0(2;) = 1 for x > 0, else 8(x ) = 0. In this article, the QCD spectral densities 
are of the type 5{ui 2 — mg), which include both the ground state and continuum state contribu¬ 
tions, I have attributed the excited state contributions to the continuum state contributions, so 
the collective continuum state contributions can be approximated as cS(lj 2 — so), then I obtain the 
result c/ (so — oj 2 ) in the hadronic representation, see Eq.(15). The doublet (D(2550), 0(2600)) 
or (Oj(2580), O}(2650)) is assigned to be the first radial excited state of the doublet (O, D*) [20!. 
The single-pole contributions come from the doublet (0(2550), 0(2600)) or (Oj(2580), O}(2650)) 

are of the form 1/ [m 2 D ^ 2 5 50 )/d( 2600 ) ~ w2 ); so the approximation c/ (so — w 2 ) is reasonable. 

Then the shifts of the masses and decay constants of the charmed-mesons can be approximated 


as 


5 m 


D/Do/D'/Di 


= 27r 


m N + m D / Do / D */ Dl 
niNm D / Do / D */ Dl 


PN<Id/D 0 /D */Di I 


( 12 ) 
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^ fD/Do - 


5fD*/D 1 = 


bpN 


D j Do m D/D 0 $ m D/D 0 


2/d/Do 7 


l D / D 0 


2m n 
( bp N 


mi 


2fD-/D 1 m D*/D 1 bm D */ Dl 


( 13 ) 


d* / Dd^R* / r 1 \2niN 

In calculations, I have used the following definitions for the decay constants of the heavy mesons 


(0|??(0)|-Do + A}) = 

(0|7, 5 (0)|£> + D) = 

<0M0)l D*+D*) = 
<0|77 5 / ,(0)| D 1 +D 1 ) = 


fPo m D 0 

m c 

f dik'd 
m c 

fD*niD*e[i , 


(14) 


q^q^ 


with summations of the polarization vectors e^A, q)e*(X, q) = —g^ + 

In the low energy limit u> —> 0, the Tr(ui, 0) is equivalent to the Born term T® orn (w, 0). Now I 
take into account the Born terms at the phenomenological side, 


T N (u 2 ) = T%°™(u 2 ) + 


with the constraint 


{ m D/D 0 /D-/D 1 u 2 ) 


+ 


rn ^D / Do / D* / D\ u2 


S 0 - W 


2 ’ 


(15) 


m 


D/D 0 /D*/D i 


TO 


D/Do/D -/Di 


— = 0 . 

SO 


(16) 


The contributions from the intermediate spin-| charmed baryon states are zero in the soft- 
limit > 0 [ 21 ] , and I only take into account the intermediate spin -1 charmed baryon states in 
calculating the Born terms, 


(D /Do/D* / Di) ( cu ) + p(uud) or n(udd) 
( D/D 0 /D*/Di) + ( cd) + p(uud) or n(udd) 


A+, £+ (cud) or YP c (cdd), 
Yi+ + (cuu) or A+, Yi+(cud) 


(17) 


where M\ c = 2.286 GeV and M^ c = 2.454 GeV |22j. I can take Mr ~ 2.4 GeV as the average 
value, where the H means either A+, E+, E+ + or E°. In the case of the bottom baryons, I take 
the approximation Mr = —./ — ~ 5.7 GeV [22] . I write down the Feynman diagrams and 

calculate the Born terms directly, and obtain the results, 

rpBonv n\ _ 2m N (m H + m N ) ( fDm 2 D g DNH 

[w 2 — (mR + TOat) 2 ] [w 2 - m 2 D ] V m c 

for the current 775(0;), 

T® orn (w, 0) -)• T® orn (w, 0) (with toat->-TO iv , D ->■ D 0 ) , ( 19 ) 

for the current p(x), 



T% OTn (co, 0) = 


for the current rjfj,(x), 


T iBorn 
N 


(w,0) 


\UJ* 


2 TOjv(mif + mjv) 

(m H + to.at) 2 ] [w 2 - m 2 D ,] 2 


(f d* tud* 9 d*nh) 


2 


—> T® orn (w, 0) (with toat —> —m at, D* —> D\) , 


( 20 ) 


( 21 ) 
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for the current 775 ^( 2 ;), where the gD/D 0 /D* /d 1 nh denote the strong coupling constants Qd/Dq/d*/d 1 na c 
and gD/D 0 /D*jD 1 NT. c - On the other hand, there are no inelastic channels for the (D/D 0 /D*/Di)°N 
and (D/Dq/D*/D i)~N interactions, and T® orn (0) = 0. In calculations, I have used the following 
definitions for the hadronic coupling constants, 


(A c /E c (p - q)\D{—q)N{p)) 
{A c /V c {p-q)\Do{-q)N(p)) 

<A C /E c (p-q)\D*(-q)N(p)) 
(A c /E c (j> - q)\D 1 (—q)N(p)) 


9a c /s c dnUa c /t. c (p - q)ilsU N (p ), 
gA c /s c D 0 NUa c /s c {p ~ q)U N {p ), 

u A c /s c (p ~ q ) ^3A C /E 0 D*JV 4 + 0-0/3 U n(p) , 

u Kc /y. c {p - q) [gA c ^ c D 1 N 4 + i ^^q^j i$u n {p ), 

( 22 ) 


where the Un and U\ a /Y, c are the Dirac spinors of the nucleon and the charmed baryons A c /E c , 
respectively. In the limit 0, the strong coupling constants ^ D , N and g p DiN have no 
contributions. 

For example, near the thresholds, the D*N can translate to the DN, D*N , 7rE c , gA c , etc, 
we can take into account the intermediate baryon-meson loops or the re-scattering effects with 
the Bethe-Salpeter equation to obtain the full D*N —> D*N scattering amplitude, and generate 
higher baryon states dynamically [23]. We can saturate the full D*N —> D*N scattering amplitude 
with the tree-level Feynman diagrams describing the exchanges of the higher resonances A c (2595), 
E c (i ), etc. While in other coupled-channels analysis, the A c (2595) emerges as a DN quasi¬ 
bound state rather than a D*N quasi-bound state [23]. The translations D*N to the ground 
states A c and E c are favored in the phase-space, as the A c (2595) and E c (i ) with J p = i have 
the average mass m_y/ « 2.7GeV 22 lM]- I n fact; m H> > s o, I can absorb the high resonances 
into the continuum states in case the high resonances do not dominate the QCD sum rules. In 
calculations, I observe that the mass-shift Stud* does not sensitive to contributions of the ground 
states A c and E c , the contributions from the spin-^ higher resonances maybe even smaller. In this 
article, I neglect the intermediate baryon-meson loops, their effects are absorbed into continuum 
contributions. 

At the low nuclear density, the condensates ( 0) PN in the nuclear matter can be approximated 
as 


(0 ] PN = + 


(23) 


based on the Fermi gas model, where the (O) and (Q)n denote the vacuum condensates and 
nuclear matter induced condensates, respectively m- I neglect the terms proportional to pp, pp, 

Pp , ■ • • at the normal nuclear matter with the saturation density = po = as the Fermi 
momentum pp = 0.27 GeV is a small quantity EH- 

I carry out the operator product expansion to the nuclear matter induced condensates (O) n 

up to dimension-5 at the large space-like region in the nuclear matter, and take into account the 
one-loop corrections to the quark condensate ( qq)N■ I insert the following term 


1 

2 ! 


ig s 


d D y^{y)^i4y)%-G a Ay) ig, 


d D z'tp(z)'y u ip( z ) — G h v { z ), 


(24) 


with the dimension D = 4 — 2e, into the correlation functions T]y(q) and T^ v {q) firstly, where the i/> 
denotes the quark fields, the G“ denotes the gluon field, the X a denotes the Gell-Mann matrix, then 
contract the quark fields with Wick theorem, and extract the quark condensate (qq)N according to 
the formula (N\q l a qp\N) = —^.{qq)Ndijbap to obtain the perturbative corrections a s (qq)N, where 
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Figure 1: The perturbative 0(a s ) corrections to the quark condensate ( qq)N■ 


the i and j are color indexes and the a and /? are Dirac spinor indexes. There are six Feynman 
diagrams make contributions, see Fig.l. Now I calculate the first diagram explicitly for the current 
i^{x) in Fig.l, 


2 T 


(c* s ,l) 


N 


(q 2 ) = - 


Tl i^-^-)(gg)Ng s g 


2 . .2e 


12 


(27r) 


D 


d°kTr i ns 4 7“ 7^ 4*75 




ft ft ' 4+ ft. — m c 


k 2 


4:Dm c (qq) N g 2 g, 2e . d (-2? ri) 


3(2 n) D dt 2n i 

Dm c (qq) N g 2 /J, 2e [1 + e(log47r - ~f E )\ 
12tt 2 


00 f d D k5 ( k 2 —t) S ((k + q) 2 — m 2 ) 
ds -^- \t-o 


s — q 


ds 


■ m n 


s — q 2 s 1 ~ e (s — m 2 ) 1+2e ’ 


(25) 


where I have used Cutkosky’s rule to obtain the QCD spectral density. There exists infrared 
divergence at the end point s = m 2 . It is difficult to carry out the integral over s, I can perform 
the Borel transform B M 2 firstly, then carry out the integral over s, 


B m2 2T, 


(«s,l) 

N 


{q 2 ) = - 


Dm c (qq) n g 2 g 2e [1 + e(log47r - j E )\ 


12t r 2 M 2 
o 

mi 


ds 


s + mi. 


s 1_e (s m 2 ) 1+2e 


exp 


(-w) 


m c (qq) n g 2 
3t r 2 M 2 6XP l M 2 


1 , „ . \ . m c (qq) N g 2 

__l 0g4T + 7E +———T 


M 2 


m c (qq) N g 2 
6t r 2 M 2 


■ exp - 


M 2 


m c (qq)Ng 2 ( m 2 \ m 2 g 2 


(26) 


where 


r°° 1 

r(0, x) = e~ x dt - e~*. (27) 

Jo t + x 

Other diagrams are calculated analogously, I regularize the divergences in D = 4 — 2e dimension, 
then remove the ultraviolet divergences through renormalization and absorb the infrared diver¬ 
gences into the quark condensate ( qq)N■ 

I calculate the contributions of other condensates at the tree level, the calculations are straight- 
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forward and cumbersome. In calculations, I use the following formulas, 


{q a {x)qp{ 0 ))jv = -t 


{qq) N + x^{qD ll q) N + ^x»x u (qD^D„q) N + 


+ (qixq)N + x^ J {q r )\D tl q) N + -x^x" {q'y\D fJ ,D„q) N H- 7 ^ 


(28) 


and 


1 Ka 


{gsqWpG^N = -^7-^ {{9sqvGq) N [(Tfjv + i{u^ v - u„ 7 M ) /t ] a/3 + {g s q /iaGq) N 

[o> i + i{u^ v - Uvhi)\ a p ~ 4 (quDuDq) N [a^ + 2i(u^ v - u v 7 M ) ^] a/3 | , 


(29) 


where 0^ = 8^- ig s ^-G a , 


{qi^q)N 

{qD^q) N 

{qj f j.D u q) N 
(qD^D^q) N 
(qixD^D^N 


= (q 'iiq) N u li , 

= ( quDq) N u M = 0 , 

= ^{qiuDq) N , 

= ^(quDuDq) N (u^u v - jgnJ] ~ \{g s qcrGq) N {u^u,, - g^ v ) , 


= 2 {qfiuDuDq) m 


1 


uyu^Uv - - (u\g^ + u^gxu + u^A/x) 


and 


-- {g s q MGq) N (uxu l _ l u„ - ma3 M ix) , 

S ab 

(GupG b ^ u )N = -g —(GG)n {ga^gpv — ga^gp^) + ^ ((® 2 + B 2 )n ) • 


(30) 


(31) 


Once analytical results at the level of quark-gluon degree’s of freedom are obtained, then I set 
ui 2 = q 2 , and take the quark-hadron duality below the continuum threshold so, and perforin the 
Borel transform with respect to the variable Q 2 = — w 2 , finally obtain the following QCD sum 
rules: 


a C a + b Cb — Cf , 

C “ = S5 eXP ('S f ) _ 5 eXP ('^)’ 


(32) 


Cb = exp 


m 2 D \ _ _so_ 
M 2 ) m 2 D 


exp 


(-&)■ 


(33) 


Cf = 


2 m N (m H + m N ) / fpm 2 D g DNH \ 


(rriH + mN ) 2 - m% \ 


M 2 


l D 


(inn + mN) 2 - mr D 


exp 


(' m H + m N f 


M 2 


- (mH + mN) 2 \ 
c {qq)N 


3 ^ M 2 

+ l< -2 11- 


log — 2r 0 


rag 

M 2 


(qU.D 0 q) N 


M 2 
4m, 


exp 


M 2 


1 - 


rag 

M 2 

mt 


2 

exp 


1 + — 
7 r 


rag 

‘ M 2 


exp 

6 - 


_m^\ 
M 2 ) 
4m 2 
3 M 2 


2 M 2 J 


( 1 qiD 0 iD 0 q) N + — { 


1 ,a s GG 


12 ' 7T 


N 


rag 

'Af 2 


exp 


(34) 



for the current 775 (re) 


Ci —> Ci (with rriN —tun , rn c —>• — m c , D —> Dq) , 


for the current rj(x), 


C a 


C b 


M 2 


exp 


M 2 




(35) 


(36) 


Cf = 


2 miv(m_y + Wat) 
(mu + tun) 2 - m 2 D , 

1 

(mn + tum ) 2 — to 


{fD^ITT'D* QD* Nh) 


D- 


1 

JP 

(' m H + m N ) 2 


m|,, - (to// + ?njv) 2 . 


exp - 


m 


D* 


M 2 


2 

+ 3 


M 2 


log — 

M 2 


2m 2 

3M 2 


r 0 


M 2 

o 

AP 


z{m) 


N 


exp 


M 2 


exp 


i + 2i 

7T 

' M 2 


4 to 2 
3 AP 


4 (qHD 0 q) N 2w 2 (g t zD 0 g)Af 2m c (qg s aGq) N 16m c (qiD 0 iD 0 q) N 


M 2 

2m 3 c (qiD 0 iD 0 q) N _ 1 a s GG 

12' 7T 


M ex p(-]v# 


3M 2 
m 2 


3M 2 


M 4 

for the current p^ix), 

Ci —> Ci (with to at —t —TOat , to, 


— to c , D* —> Di) , 


(37) 


(38) 


for the current 775^(2;), where * = a,b,f. In this article, I neglect the contributions from the 
heavy quark condensates ( QQ ), (QQ) = — 12 ^ mQ { aaGG ) up to the order 0(a s ) (here I count the 

condensate ( a “ GG ) as of the order 0(a s )), the heavy quark condensates have practically no effect 
on the polarization functions, for detailed discussions about this subject, one can consult Ref. [9]. 
In Ref. [25], Buchheim, Hilger and Kampfer study the contributions of the condensates involve the 
heavy quarks in details, the results indicate that those condensates are either suppressed by the 
heavy quark mass toq or by the additional factor (or g 2 /( 47r) 2 ). Neglecting the in-medium 
effects on the heavy quark condensates cannot affect the predictions remarkably, as the main 
contributions come from the terms ( qq)N■ 

Differentiate above equation with respect to r = then eliminate the parameter b (a), I can 
obtain the QCD sum rules for the parameter a (&), 

Cf(-&)C b -C b (-£)C f 
Ca(-£)C b -C b {-£)C a ’ 

Cf(-&)C a -C a (-£)C f 

c b {-£)c a -c a {-l)c h - 

With the simple replacements m c —t m b , D/D 0 /D*/D\ —> B/B 0 /B*/Bi , A c —> A b and £ c —> T, b , 
I can obtain the corresponding the QCD sum rules for the bottom mesons in the nuclear matter. 


a = 

b = 


3 Numerical results and discussions 

9 3 

At the normal nuclear matter with the saturation density pw = Po = gyr, where the Fermi 
momentum pp = 0.27 GeV is a small quantity, the condensates ( 0) PN in the nuclear matter can 
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be approximated as ( 0) PN = (O) + ^~{G)n, the terms proportional to p F , p F , p F , ■ ■ ■ can be 
neglected safely, where the (G) = (O|0|O) and (O)n = {N\G\N) denote the vacuum condensates 
and nuclear matter induced condensates, respectively HU- 

The input parameters at the QCD side are taken as pn = (0.11 GeV) 3 , (qq)N = — (2m n), 

r° G ) N = — 0.65 GeV(2mjv), &n = 45MeV, m u + rrid = 12MeV, ( qUD 0 q)N = 0.18 GeV(2mjv), 
( qg s aGq)N = 3.0GeV 2 (2mjv), (qiD 0 iD 0 q)pf + |( qg s aGq)N = 0.3 GeV 2 (2mjv), tun = 0.94GeV 
I TU] . m c = (1.3 ± 0.1) GeV, nil, = (4.7 ± 0.1) GeV, a s = 0.45 and p = 1 GeV. If we take the 
normalization {N(p)\N(p')) = (2 tt) 3 5 3 (p — p'), then ( 0) PN = (O) + Pn(G)n, the unit 2 toat 
in the brackets in the values of the condensates ( qq)N , ( a ‘ GG )N, disappears. I choose the 
values of the nuclear matter induced condensates determined in Ref.|10j. which are still widely 
used in the literatures. Although the values of some condensates are updated, those condensates 
are irrelevant to the present work. The updates focus on the four-quark condensate (261 . In this 
article, I take into account the condensates up to dimension-5, the four-quark condensates have 
no contributions, the dominant contributions come from the nuclear matter induced condensate 
(qq) N , {qq)N = — a + md (2mjv)- The value m u + m<i = 12MeV is obtained from the famous Gell- 
Mann-Oakes-Renner relation at the energy scale p = 1 GeV, while the value un = 45 MeV is still 
widely used [25] , 

The parameters at the hadronic side are taken as m B = 1.870 GeV, m B = 5.280 GeV, m Bo = 
2.355 GeV, m Bo = 5.740 GeV, m D . = 2.010 GeV, m B > = 5.325 GeV, m D] = 2.420 GeV, m Bl = 
5.750 GeV, f D = 0.210 GeV, f B = 0.190 GeV, f Do = 0.334^ GeV, f Bo = 0.280^ GeV, f D . = 


0.270 GeV, f B , = 0.195 GeV, f Dl = 0.305 GeV, f Bl = 0.255 GeV, s° D = (6.2 ± 0.5) GeV 2 , s° B = 
(33.5 ± 1.0) GeV 2 , s° D , = (6.5 ± 0.5) GeV 2 , s° B , = (35.0 ± 1.0) GeV 2 , s° Do = (8.0 ± 0.5) GeV 2 , 
4 0 = (39.Oil.0) GeV 2 , s° Di = (8.5±0.5) GeV 2 and s° Bi = (39.Oil.0) GeV 2 , which are determined 
by the conventional two-point correlation functions using the QCD sum rules mm - 1 neglect 
the uncertainties of the decay constants to avoid double counting as the main uncertainties of the 
decay constants originate from the uncertainties of the continuum threshold parameters Sq. 

The value of the strong coupling constant gDNA c is gK c DN = 6.74 from the QCD sum rules 
58 ] , while the average value of the strong coupling constants gA c DN and gs c DN from the light-cone 
QCD sum rules is 9A = ol *+ RD11 = 6.775 [29], those values are consistent with each other. The 
average value of the strong coupling constants gA c D*N and gs c D*N from the light-cone QCD sum 
rules is 9A ^ DtN = 3.86 [29]. In this article, I take the approximation gDNA c ~ gDNH c ~ 

gBNA b ~ gBNY, b ~ gD 0 NA c ~ gD 0 NS c ~ gB 0 NA b ~ gB 0 NS b ~ 6.74 and gA c D*N ~ gY. c D*N ~ 
gA c D b N ~ gn c D 1 N ~ gA b B*N ~ g?. b B*N ~ gA b B ± N ~ g^2 b B 1 N ~ 3.86. 

In Figs.2-3, I plot the shifts of the masses and decay constants of the heavy mesons in the 
nuclear matter with variations of the Borel parameter M 2 , respectively. From the figures, I can 
see that there appear platforms. In this article, I choose the Borel parameters M 2 according to 
the criterion that the uncertainties originate from the Borel parameters M 2 are negligible. The 
values of the Borel parameters M 2 are shown explicitly in Table 1. From Figs.2-3 and Table 1, 
I can see that the Borel parameters M 2 in the QCD sum rules for the mass-shift 5m and decay- 
constant-shift 5f of the same meson are different. It is not un-acceptable, as the mass-shift 6m 
and decay-constant-shift 5f come from different QCD sum rules, not coupled QCD sum rules, see 
Eq.(39), the platforms maybe appear in different places in different QCD sum rules. 

I can obtain the shifts of the masses and decay constants of the heavy mesons in the nuclear 
matter in the Borel windows, which are shown explicitly in Table 2. From the Table 2, I can 
obtain the fractions of the shifts < 5%, < 10% Sm B/B , /B()/Bl = 

m D/D , /Do/Dj - JD/D-/D 0 /D 1 - rn B / B* / Bq / Bl 


(5 — 15)% and ‘ > ~f B/B * /B ° /£il = (25 — 55)%, which are shown explicitly in Table 3. In calculations, 

I observe that the main contributions come from the terms m c (qq)]y and mb(qq)N- From Table 
3, I can see that the next-to-leading order corrections a s (qq)N are important. In the case of 
the shifts <5 " 1J3 * /Bl and ; the next-to-leading order contributions a s {qq)N and the leading 

' ^B /B\ JB /B\ 


order contributions (qq)N are almost equivalent. In this article, I choose the special energy scale 
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Figure 2: (Color online) The shifts of the masses of the heavy mesons in the nuclear matter with 
variations of the Borel parameter M 2 , the I (II) denotes contributions up to the next-to-leading 
order (leading order) are included. 
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Figure 3: (Color online) The shifts of the decay constants of the heavy mesons in the nuclear 
matter with variations of the Borel parameter M 2 , the I (II) denotes contributions up to the 
next-to-leading order (leading order) are included. 
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Stud 

5m,D * 

5niD 0 

SrriD i 

Sms 

6ms * 

SrriBo 

Stub x 


4.4 -5.4 

4.6-5.6 

6.0- 7.0 

6.6-7.6 

29-33 

30-34 

32-36 

32-36 


SfD 

5/d* 

5/d 0 

Sf Dl 

Sf B 

SfB * 

S/b 0 

Sf Bl 

~M^ 

1.9- 2.9 

3.5-4.5 

4.3-5.3 

5.3-6.3 

25-29 

27-31 

30-34 

31-35 


Table 1: The Borel parameters in the QCD sum rules for the shifts of the masses and decay 
constants of the heavy mesons in the nuclear matter, the unit is GeV 2 . 



SniD 

Stud* 

SniD 0 

SmDx 

Sms 

Stub * 

SmB 0 

SmBi 

NLO 

-72 

-102 

80 

97 

-473 

-687 

295 

522 

LO 

-47 

-70 

54 

66 

-329 

-340 

209 

260 

m 

-48 








m 

+45 




+60 




HD 

-46 




-242 





SfD 

5/d* 

5/d 0 

Sf Dl 

SfB 

5/b* 

5/b 0 

SfB i 

NLO 

-6 

-26 

11 

31 

-71 

-111 

56 

134 

LO 

-4 

-18 

7 

21 

-48 

-55 

39 

67 

na 

-2 




-23 





Table 2: The shifts of the masses and decay constants of the heavy mesons in the nuclear matter, 
where the NLO (LO) denotes contributions up to the next-to-leading order (leading order) are 
included, the unit is MeV. 


/r = 1 GeV. The logarithm log ^ in the next-to-leading contributions is very large and enhances 

the next-to-leading contributions greatly. Although the nuclear matter induced condensates evolve 

with the renormalization group equation, their evolving behaviors with the energy scales are not 

well known, as this subject has not been studied in details yet at the present time. A larger energy 

2 

scale /i can lead to smaller logarithm log T -jff therefore more reasonable predictions. In Table 4, I 
present the main uncertainties, which originate from the uncertainties of the heavy quark masses 
and the continuum threshold parameters. 

The mass-shifts of the negative (positive) parity mesons are negative (positive), the decays 
of the high charmonium states to the DD and D*D* (DqDq and DiD\) pairs are enhanced 
(suppressed) in the phase space, and we should take into account those effects carefully in studying 
the production of the J/ip so as to identifying the quark-gluon plasmas. 



8m d 

8m £>* 

4mo 0 

dmD 1 

8m b 

8m b* 

SmB 0 

SmB 1 


mD 

m D * 

mo n 

mDi 

ms 

m B * 

mB 0 

msi 

NLO 

-4% 

-5% 

3% 

4% 

-9% 

-13% 

5% 

9% 

LO 

-3% 

-3% 

2% 

3% 

-6% 

-6% 

4% 

5% 


Sf D 

S.fn* 

SfD o 

SfD i 

SfB 

Sfn* 

SfB 0 

SfB 1 


Id 

Id* 

fa n 

i 


/b* 

/s n 

f* 1 

NLO 

-3% 

-10% 

6% 

10% 

-37% 

-57% 

24% 

53% 

LO 

-2% 

-7% 

4% 

7% 

-25% 

-28% 

17% 

26% 


Table 3: The fractions of the shifts of the masses and decay constants of the heavy mesons in 
the nuclear matter, where the NLO (LO) denotes contributions up to the next-to-leading order 
(leading order) are included. 
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5{5mp ) 

5{5mp*) 

5(6mp 0 ) 

5(5mp 1 ) 

5 (5mp) 

5(5mp *) 

5(5mp 0 ) 

5(5mp 1 ) 

5itiq 

±14 

±4 

±26 

±6 

±18 

±1 

±25 

±1 

5sq 

±9 

±14 

±12 

±13 

±65 

±80 

±39 

±70 


5{5fp) 

5(5 fp.) 

8(5 f do) 

8(5 f Dl ) 

8(8 f B ) 

5(5 fpQ 

5(8 fp 0 ) 

5(5 f Bl ) 

5m Q 

±1 

±1 

±3 

±2 

±2 

±1 

±3 

±0 

5s 0 

±1 

±7 

±4 

±8 

±21 

±25 

±15 

±34 


Table 4: The uncertainties of the shifts of the masses and decay constants of the heavy mesons 
in the nuclear matter originate from the uncertainties of the heavy quark masses and continuum 
threshold parameters, where the unit is MeV. 



ap 

ap * 

ap o 

ap i 

ap 

a B * 

&Bq 

ap i 

NLO 

-1.1 

-1.5 

1.3 

1.6 

-8.9 

-12.9 

5.6 

9.9 

LO 

-0.7 

-1.1 

0.9 

1.1 

-6.2 

-6.4 

4.0 

5.0 


Table 5: The heavy-meson-nucleon scattering lengths, where the NLO (LO) denotes contributions 
up to the next-to-leading order (the leading order) are included, the unit is fm. 


The currents Qq and < 3*759 (also Qj p q and Q7 m 759) are mixed with each other under the 
chiral transformation q —> e lai5 q, the currents Qq , 0*759, Q^^q, Qj^^q are not conserved in 
the limit m q —> 0, it is better to take the doublets ( D,D 0 ) and (D*,D i) as the parity-doublets 
rather than the chiral-doublets. The quark condensate {qq) pN serves as the order parameter, and 
undergoes reduction in the nuclear matter, the chiral symmetry is partially restored; however, 
there appear new medium-induced condensates, which also break the chiral symmetry. In this 
article, the (O)n are companied by the heavy quark masses toq, uiq or m,Q, the net effects cannot 
warrant that the chiral symmetry is monotonously restored with the increase of the pjv. When 
the pat is large enough, the order parameter (qq) PN 0, the chiral symmetry is restored, the 
Fermi gas approximation for the nuclear matter breaks down, and the parity-doublets maybe have 
degenerated masses approximately. In this article, I study the parity-doublets at the low pjv, the 
mass breaking effects of the parity-doublets maybe even larger, see Table 2. We expect that smaller 
mass splitting of the parity-doublets at the high nuclear density is favored, however, larger mass 
splitting of the parity-doublets at the lower nuclear density cannot be excluded. In Refs. mm 
the mass center Trip of the pseudoscalar mesons increases in the nuclear matter while the mass 
center mg of the scalar mesons decreases in the nuclear matter, the mass breaking effect ms — Trip 
of the parity-doublets is smaller than that in the vacuum. 

In Table 5, I show the scattering lengths ap, ap*, ap 0 , ap 1 , ap, apt, «b 0 i a Bl explicitly, the 
ap, ap *, ap, ap » are negative, which indicate the interactions DN, D*N, BN, B*N are attractive, 
the ap 0 , flfli, ap 0 , ap 1 are positive, which indicate the interactions DqN, D\N , BqN , BqN are 
repulsive. It is difficult (possible) to form the DqN , D\N, BqN, BqN {DN, D*N, BN, B*N) 
bound states. In Ref. m, Zhang studies the S'-wave DN and D*N bound states with the QCD 
sum rules, the numerical results indicate that the £ c (2800) and A c (2940) can be assigned to be 
the .S'-wave DN state with J p = ~ and the S'-wave D*N state with J p = § , respectively. 

In the present work and Refs. iii, the correlation functions are divided into a vacuum 
part and a static one-nucleon part, and the nuclear matter induced effects are extracted explicitly; 
while in Refs. (muz], the pole terms of (or ground state contributions to) the hadronic spectral 
densities of the whole correlation functions are parameterized as All (oj ) = F + 5(uj — m+) — F_ 5 (w ± 
m_), where m± = m ± Am and F± = F ± A F, and QCD sum rules for the mass center rn and 
the mass splitting Am are obtained. In the leading order approximation, the present predictions 
of the 5mp and 5mp are compatible with that of Refs. m fl5] and differ greatly from that of 
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Figure 4: (Color online) The contributions of the perturbative term (I) and quark condensate term 
(II) in the QCD sum rules for the D mesons in the vacuum. Furthermore, I show the mass mu (III) 
and decay constant fo (IV) explicitly, which are normalized to be 1 at the value M 2 = 1.5GeV 2 . 


Refs, pa QI], see Table 2. The values obtained from the QCD sum rules depend heavily on the 
Borel windows, the values extracted from different Borel windows especially in different QCD sum 
rules maybe differ from each other greatly. 

In Refs. [161117) . the authors study the masses of the heavy mesons in the nuclear matter directly 
by including both the vacuum part and the static one-nucleon part in the QCD sum rules, then 
the continuum contributions are well approximated by Pqcd(u 2 )9 (lu 2 — w^ 2 ^, where the u>± 2 are 
the continuum threshold parameters, it is one of the advantages of Refs. PHUT!. However, they 
define the moments S n (M 2 ) to study the mass-shifts, 


S n {M 2 ) 


dujoj n An(w) exp 



(40) 


the odd moment o = Sq(M 2 ) and the even moment e = S\{M 2 ), then obtain = —S 2 (M 2 ) 

and jTnjjTj = —S'i(M 2 ) by assuming the F± and m± are independent on the Borel parameters at 
the phenomenological side. In fact, d ^/° M ‘i) ^ — S 2 (M 2 ) and jrnjfTj —.^(M 2 ) at the operator 
product expansion side according to the QCD spectral densities An(w), which depend on the Borel 
parameters explicitly, the approximations = —S 2 (M 2 ) and jnjj jt; = ~S 3 (M 2 ) lead to 

undetermined uncertainties. 

In Refs. [16] [T7j, the perturbative 0(a s ) corrections to the perturbative terms are taken into 
account. In the QCD sum rules for the pseudoscalar D mesons in the vacuum, if we take into 
account the perturbative 0(a s ) corrections to the perturbative term and vacuum condensate term, 
the two criteria (pole dominance and convergence of the operator product expansion) of the QCD 
sum rules leads to the Borel window M 2 = (1.2 — 1.8) GeV 2 , the resulting predictions of the 
mass mo and decay constant fu are consistent with the experimental data. In Fig.4, I plot 
the contributions of the perturbative term and quark condensate term in the operator product 
expansion. From the figure, I can see that the main contributions come from the perturbative 
term, the quark condensate ( qq ) plays a less important role. 

The modifications of the condensates in the nuclear matter are mild, for example, (qq) pN ~ 
0.64 (qq), while the perturbative contributions are not modified (or modified slightly by introduc¬ 
ing a minor splitting Aso, = so ± Aso) by the nuclear matter. If we turn on the in-medium 
effects, the contributions of the quark condensate are even smaller, the Borel windows are deter- 
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mined dominantly by the perturbative terms mm- If the perturbative 0(a 2 ) corrections to the 
perturbative terms are also included, the contributions of the perturbative are even larger EH, the 
QCD sum rules are dominated by the perturbative terms, which are not (or slightly) affected by 
the nuclear matter. It is not favored to extract the mass-shifts in the nuclear matter, and impairs 
the predictive ability. 

In the present work and Refs. min ng, the correlation functions are divided into the vacuum 
part and the static one-nucleon part, which are of the orders 0(0) and 0(pn ), respectively. We 
can obtain independent QCD sum rules from the two parts respectively. The QCD sum rules 
correspond to the orders 0(0) and 0{pn) respectively can have quite different Borel parameters. 
In this article, I separate the nuclear matter induced effects unambiguously, study the QCD sum 
rules correspond to the order 0(pjv), and determine the Borel parameters by the criteria of the 
QCD sum rules. 

In the conventional QCD sum rules, we usually choose the Borel parameters M 2 to satisfy the 
following three criteria: 

1. Pole dominance at the phenomenological side; 

2. Convergence of the operator product expansion; 

3. Appearance of the Borel platforms. 

In the present work and Refs. mm EE the nuclear matter induced effects are extracted 
explicitly, the resulting QCD sum rules are not contaminated by the contributions of the vacuum 
part, the Borel windows are determined completely by the nuclear matter induced effects, it is 
the advantage. As the QCD spectral densities are of the form S(uj 2 — TOq), we have to take the 
hadronic spectral densities to be the form 5(lu 2 — m 2 H ) and model the continuum contributions 
with the function <5(w 2 — So), and determine the s o by some constraints, see Eq.(16), where the H 
denotes the ground state and excited state heavy mesons. In this article, I attribute the higher 
excited states to the continuum contributions, the (5-type hadronic spectral densities make sense. So 
the pole dominance at the phenomenological side can be released as the continuum contributions 
are already taken into account. Furthermore, I expect that the couplings of the interpolating 
currents to the excited states are more weak than that to the ground states, the uncertainties 
originate from continuum contributions are very small. For example, the decay constants of the 
pseudoscalar mesons 7t(140) and 7r(1300) have the hierarchy /^( 1300 ) <C /,r(i 40 ) from the Dyson- 
Schwinger equation [30,, the lattice QCD [31], the QCD sum rules [32], etc, or from the experimental 
data [33] . 

In the present work and Refs. [Hill 115] , large Borel parameters are chosen to warrant the 
convergence of the operator product expansion and to obtain the Borel platforms, and small Borel 
parameters cannot lead to platforms. In the Borel windows, where the platforms appear, the main 
contributions come from terms ( qq)N , the operator product expansion is well convergent. The 
criteria 2 and 3 can be satisfied. The continuum contributions are not suppressed efficiently for 
large Borel parameters compared to that for small Borel parameters. In calculations, I observe 
that the predictions are insensitive to the so, the uncertainties originate from the continuum 
threshold parameters So are very small in almost all cases, the large Borel parameters make sense. 
Furthermore, the continuum contributions are already taken into account. On the other hand, from 
Eqs.(8-9) and Eqs.(12-13), we can see that the mass-shifts Suid/Dq/d-/d x and decay constant shifts 
<5/ d/Do/d */Di reduce to zero in the limit p^ —>• 0, the QCD sum rules correspond to the nuclear 
matter induced effects decouple, their Borel parameters (irrespective of large or small) are also 
irrelevant to the ones in the QCD sum rules for the vacuum part of the correlation functions. So 
the present predictions are sensible. 

The predictions depend on the in-medium hadronic spectral functions [34] . for example, there 
are two generic prototypes of the in-medium spectral functions for the p meson, they differ in 
details at the low mass end of the spectrum. The Klingl-Kaise-Weise spectral function emphasizes 
the role of chiral in-medium 7T7t interactions [35] , while the Rapp-Wambacli spectral function 
focuses on the role of nucleon-hole, A(1232)-hole and iV*(1520)-hole excitations [36]. Both of the 
spectral functions account quite well for the low-mass enhancements observed in dilepton spectra 
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from high-energy nuclear collisions. However, the QCD sum rules analysis of the lowest spectral 
moments reveals qualitative differences with respect to their Brown-Rho scaling properties |34j . 
If the simple spectral densities FS( w 2 — M|y y ) analogous to the ones in Refs.[Thl 17] are taken, 
where the P denotes the pseudoscalar mesons 7r, p c , the V denotes the vector mesons p, w, <f>, J/tjj, 
the F denotes the constant pole residues, the in-medium mass-shifts 5M P / V are smaller than zero 
qualitatively 122 • i expect that the S-wave mesons q'q , cq, cc with the spin-parity J p = 0 (or 
1 _ ) have analogous in-medium mass-shifts, at least qualitatively. Further studies based on more 
sophisticated hadronic spectral densities are needed. 

In fact, there are controversies about the mass-shifts of the D and B mesons in the nuclear 
matter, some theoretical approaches indicate negative mass-shifts [23], while others indicate pos¬ 
itive mass-shifts [38] . The different predictions originate mainly from whether or not the heavy 
pseudoscalar and heavy vector mesons are treated on equal footing in the coupled-channel ap¬ 
proaches. If we obtain the meson-baryon interaction kernel by treating the heavy pseudoscalar 
and heavy vector mesons on equal footing as required by heavy quark symmetry, the mass-shift 
SM d is negative [23], which is consistent with the present work; furthermore, the attractive D- 
nucleus interaction can lead to the formation of D-nucleus bound states, which can be confronted 
to the experimental data in the future directly [351 . 

The upcoming FAIR project at GSI provides the opportunity to study the in-medium properties 
of the charmoniums or charmed hadrons for the first time, however, the high mass of charmed 
hadrons requires a high momentum in the antiproton beam to produce them, the conditions for 
observing in-medium effects seem unfavorable, as the hadrons sensitive to the in-medium effects 
are either at rest or have a small momentum relative to the nuclear medium. We have to find 
processes that would slow down the charmed hadrons inside the nuclear matter, but this requires 
more detailed theoretical studies. Further theoretical studies on the reaction dynamics and on the 
exploration of the experimental ability to identify more complicated processes are still needed. 

4 Conclusion 

In this article, I divide the two-point correlation functions of the scalar, pseudoscalar, vector and 
axialvector currents in the nuclear matter into two parts, i.e. the vacuum part and the static 
one-nucleon part, then study the in-nredium modifications of the masses and decay constants by 
deriving QCD sum rules from the static one-nucleon part of the two-point correlation functions. 
In the operator product expansion, I calculate the contributions of the nuclear matter induced 
condensates up to dimension 5, especially I calculate the next-to-leading order contributions of 
the in-nredium quark condensate and obtain concise expressions, which also have applications in 
studying the mesons properties in the vacuum. In calculation, I observe that the next-to-leading 
order contributions of the in-medium quark condensate are very large and should be taken into 
account. 

All in all, I study the properties of the scalar, pseudoscalar, vector and axialvector heavy 
mesons with the QCD sum rules in a systematic way, and obtain the shifts of the masses and 
decay constants in the nuclear matter. The numerical results indicate that the mass-shifts of the 
negative parity and positive parity heavy mesons are negative and positive, respectively. For the 
pseudoscalar meson D, I obtain the prediction SMp < 0, which is in contrast to the prediction 
in Refs. mm, where the mass-shift is positive 5M P > 0. In Refs. [HE], the authors study 
the masses of the heavy mesons in the nuclear matter directly by including both the vacuum part 
and static one-nucleon part with the QCD sum rules, and parameterize the spectral density of the 
whole correlation functions by a simple function An(w) = F + S(to — m+) — F-5(u> + m_). I discuss 
the differences between the QCD sum rules in the present work and that in Refs. ng hi] in details, 
and show why I prefer the present predictions. In the present work and Refs. [HI III E2 IHl E|, 
the finite widths of the mesons in the nuclear matter are neglected, further studies based on the 
more sophisticated hadronic spectral densities by including the finite widths are needed. 
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As the masses of the heavy meson paries, such as the DD , D*D *, DqDq, D\D\ are modified 
in the nuclear environment, we should take into account those effects carefully in studying the 
production of the J/'ip (and T) so as to identifying the quark-gluon plasmas. Furthermore, I study 
the heavy-meson-nucleon scattering lengths as a byproduct, and obtain the conclusion qualitatively 
about the possible existence of the heavy-meson-nucleon bound states. 
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